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ABSTRACT 

We investigate a model of disk galaxies whereby viscous evolution of the gaseous 
disk drives material inwards to form a proto-bulge. We start from the standard picture 
of disk formation through the settling of gas into a dark halo potential well, with the 
disk initially coming into centrifugal equilibrium with detailed conservation of angular 
momentum. We derive generic analytic solutions for the disk-halo system after adia- 
batic compression of the dark halo, with free choice of the input virialized dark halo 
density profile and of the specific angular momentum distribution. We derive limits on 
the final density profile of the halo in the central regions. Subsequent viscous evolution 
of the disk is modelled by a variation of the specific angular momentum distribution 
of the disk, providing analytic solutions to the final disk structure. The assumption 
that the viscous evolution timescale and star formation timescale are similar leads 
to predictions of the properties of the stellar components. Focusing on small 'expo- 
nential' bulges, ones that may be formed through a disk instability, we investigate 
the relationship between the assumed initial conditions, such as halo 'formation', or 
assembly, redshift z/, spin parameter A, baryonic fraction F, and final disk properties 
such as global star formation timescale, gas fraction, and bulge-to-disk ratio. We find 
that the present properties of disks, such as the scale length, are compatible with a 
higher initial formation redshift if the re-distribution by viscous evolution is included 
than if it is ignored. We also quantify the dependence of final disk properties on the 
ratio F/X, thus including the possibility that the baryonic fraction varies from galaxy 
to galaxy, as perhaps may be inferred from the observations. 

Key words: galaxies: formation — galaxies: structure — galaxies: spiral — cosmol- 
ogy: theory — dark matter 



1 INTRODUCTION 

The current picture of disk galaxy formation and evolution 
has as its basis the dissipative infall of baryons within a 
dominant dark halo potential well (White & Rees 1978). 
The collapse and spin-up of the baryons, with angular mo- 
mentum conservation, can provide an explanation for many 
of the observed properties of disks, with the standard ini- 
tial conditions of baryonic mass fraction _F ~ 0.1 and dark 
halo angular momentum parameter A ~ 0.07 (Fall & Efs- 
tathiou 1980; Gunn 1982; Jones & Wyse 1983; Dalcanton, 
Spergel & Summers 1997; Hernandez & Gilmore 1998; Mo, 
Mao & White 1998; van den Bosch 1998). Galaxies such as 
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the Milky Way which have an old stellar population in the 
disk, must, within the context of a hierarchical-clustering 
scenario, evolve through only quiescent merging/accretion, 
so as to avoid excessive heating and disruption of the disk 
(Ostriker 1990). Further, the merging processes with signif- 
icant substructure cause angular momentum transport to 
the outer regions, which must somehow be suppressed to al- 
low the formation of extended disks as observed (e.g. Zurek, 
Quinn & Salmon 1988; Silk & Wyse 1993; Navarro & Stein- 
metz 1997). 

Thus here we adopt the simplified picture that disk 
galaxies form from smooth gaseous collapse to centrifugal 
equilibrium, within a steady dark halo potential. We discuss 
where appropriate below how this may be modified to take 
account of subsequent infall, or earlier star formation. Our 
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model incorporates the adiabatic response of the dark halo 
to the disk infall, and we provide new, more general, analytic 
solutions for the density profile, given a wide range of initial 
density profiles and angular momentum distributions. We 
provide new insight into the 'disk-halo' conspiracy within 
the context of this model, demonstrating how an imperfect 
conspiracy is improved by the disk-halo interaction. We ex- 
plicitly include subsequent viscous evolution of the gas disk 
to provide the exponential profile of the stellar disk, and 
develop analytic expressions that illustrate the process. The 
resulting radial inflow builds up the central regions of the 
disk and we investigate the properties of 'bulges' that may 
form as a consequence of instabilities of the central disk. 
We derive new constraints on the characteristic redshift of 
disk star formation. We obtain a simple relation connecting 
the initial conditions, such as spin parameter and baryonic 
mass fraction, to the efficiency of viscous evolution and star 
formation. 



2 THE DISK GALAXY FORMATION MODEL 

In this section we shall derive the mass profiles of disk and 
halo after the collapse of the baryons. We shall follow earlier 
treatments of disk galaxy formation (e.g. Mo, Mao & White 
1998) by assuming that the virialized dark halo, mixed with 
baryonic gas, is 'formed' - or at least assembled - at red- 
shift Zf. This virialized halo has a limiting radius r2oo within 
which the mean density is 200pcrit{zf), and contains a bary- 
onic mass fraction F. Then 



^200 = 



V20C 



Mtot = 



G 



WGH{zf)' 



(1) 



where H{zf) is the value of the Hubble parameter at redshift 
Zf, Mtot{zf) is the total mass within virialized radius r2oo, 
and V200 is the circular velocity at r2oo- 

The baryonic gas cools and settles into a disk, caus- 
ing the dark halo to contract adiabatically (Blumenthal et 
al. 1986). The specific angular momentum distribution of 
the gas is assumed to be conserved during these stages. We 
shall include below the subsequent re-arrangement of the 
disk due to angular momentum transport. This we investi- 
gate by variation of the disk angular momentum distribution 
function, choosing an appropriate analytic functional form. 
Let md{r) and mh{r) respectively denote the fraction of the 
total baryonic mass, and total dark mass, that is contained 
within radius r, and denote the baryonic mass angular mo- 
mentum distribution function by 



"^d (< j) = f U/jmax) , 



(2) 



where jmax is the maximum specific angular momentum of 
the disk. 

We will be requiring that the functional form, 

f {3/ jmax), vary as the disk evolves, and it is convenient 
to introduce the notation i = j/jmax and define 



Jo 



(3) 



which represents the area above the angular momentum dis- 
tribution function curve /(£) for < < 1. We will mimic 
the effects of viscous evolution by decreasing the value of c/ 
in our evolving disk models in section 4 below. 

In terms of this parameter the total disk angular mo- 
mentum is: 

Jo 



Mdjr, 



f{£)dl] = FMt (4) 



with Md = FMtot- Thus jmaxCf is the average specific an- 
gular momentum of the disk material. 

The specific angular momentum of the disk material 
is assumed to follow that of the dark halo, but in gen- 
eral will not be a simple analytic function (e.g. Quinn & 
Binney 1992). For illustration, we adopt an analytic mono- 
tonic increasing function f{b,i) containing a free parameter 
b, with < b < 1. We require the initial angular momentum 
distribution to be scale free, representing the angular mo- 
mentum distribution of the virialized halo, and will adopt 
f{b — 0,£) = We shall vary the value of the parameter 
b to mimic the effects of viscous evolution on the angular 
momentum distribution. 

The total energy is: 



2r2oo 



eoMtotV2oo 



(5) 



where eo is a constant of order unity, depending on the dark 
halo density profile, and since it is constant for any specific 
halo model, we can take eo = 1 without loss of generality. 
The spin parameter A is by definition 



X = Jt. 



\'/'G-'M-r. 



(6) 



Thus the mean specific angular momentum of the disk ma^ 
terial may be expressed as 



(7) 



Assuming spherical symmetry, the rotationally-supported 
disk has mass profile given by 



. , I \/GMtot{md)r{md) \ 
md = f {j/jmax) = f -. = fW, (8) 



while the initial virialized halo mass profile is 

g{Rini) = mh{Rini) = Mini(Ri„i)/Mtot, (9) 

with R = r/r2oo. 



2.1 Constraints on the Final Dark Halo Profile 
and Mass Angular Momentum Distribution 
Function 

The above equations describe the disk and halo just upon 
the settling of the gas disk to the mid-plane, prior to the sub- 
se(jueiit adiabatic compression of the halo. A self-consistent 
calculation of the modified disk and halo density profiles 
may be made by consideration of the adiabatic invariance of 
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the angular action, 1$ = J ve ■ rdd = \J GMtot{'r)r, together 
with the assumption of no shell crossing (cf. Blumenthal et 
al. 1986). 

Suppose a dark matter particle initially at rini finally 
settles at r{md), the radius within which the dark halo mass 
fraction is nih- Then under adiabatic invariance the disk 
mass profile, rrid, and the halo mass profile, ruh, are related 
by: 



GMtot (md) r (md) 



GMtot{mh)r{mh) (10) 
= GMini (nih) Tini {rrih) (11) 
= GMtotmhr2Qog^^ {rrih) , (12) 

where g^^{rnh) is the inverse function of g{R), the initial 
virializcd halo mass profile. 

Further manipulation of these relations is simplified by 
introduction of the parameter ^, given by 



\/GMtotr2m _ Cf 



Jmax 



(13) 



Generally ^ is a quantity that is closely related to the overall 
disk collapse factor. 

Prom equations (8) and (12), we have 



1/2 



(14) 
(15) 



where I is the normalized specific angular momentum. 
Again, < i' < 1, and t = \ corresponds to the maximum 
specific angular momentum of the disk, which occurs at the 
edge of the disk, equivalently at the disk cutoff radius. The 
fraction of the dark matter contained within the disk thus 
has a maximum value, mhc, given by ^ = 1 in the above 
equation, and for radii with mh > Tn/ic, nid = 1. 

To illustrate the physical meaning of these parame- 
ters, consider the rigid singular isothermal halo, for which 
girrihc) = TUhc = Rc- Then from equation (15), Rc = R{£ = 
1) — 1/^ which corresponds to the disk cutoff radius. Thus 
in this case, remembering that R is the normalized radius, 
^ = l/_Rc is the disk collapse factor. 

Up to now we know the mass profile of disk and halo 
after collapse, in terms of the normalized specific angular 
momentum £, as given in equations (14) and (15), for given 
forms of the angular momentum distribution function, /, 
and initial virialized dark halo mass profile, g. Next we shall 
obtain the relation between i and radius R, to complete the 
derivation of the mass profiles of disk and halo after collapse. 

Returning to a general halo density profile, the total 
mass contained within the radius corresponding to rUh is 



M(mh) = Mtot ((1 - Fynih + Fnid) . 
From equations (9) - (12) and (16), we have 
GMini{mh)rini{mh) 



R = 



GM{mh)r2oo 
mhg~^{mh) 



(16) 

(17) 
(18) 



{1 - F)mh + Frud' 
Introducing the radius variable x and coefficient co as 

X = ^(1-m (19) 
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CO = ^F/{l-F), (20) 

we may finally derive the functional dependences on £ of the 
disk mass rud, of the halo mass ruh, and of the radius x: 



rud 



(21) 
(22) 



(23) 



^mh{e) + Co f{i)' 

Thus £ can be thought of as a normalized radius. As we 
shall see later, co is a measure of the compactness of the 
final collapsed disk due to the competition between the spin 
parameter A and the baryonic mass fraction F. 

These equations (21) - (23) can be used to derive disk 
and halo properties for a free choice of virialized halo pro- 
file g{R) and angular momentum distribution function f{i). 
Within the disk cutoff radius, with < £ < 1, the disk 
surface density, circular velocity and the disk-to-dark mass 
ratio as function of radius £ ot R have the generic forms: 

10H{z)FV2oo 1 df d£ 
2-kG Rd£dR' 



^R ' 
cofi£) 



(24) 
(25) 



where Mdi£) = Md{£ = l)md(,£) = Mdmd{£) and Mh{£) is 
defined similarly. 

The circular velocity at radii beyond the disk cutoff, 
but within the halo, is given by: 



V200 \/mhg~^{mh) 
R ' 
mhg~^{mh) 



(27) 

^ - il-F)m, + F- ^^^> 

Armed with these relations, one may now look at var- 
ious initial virialized lialo density profiles and angular mo- 
mentum distributions, and determine the allowed parameter 
space from observed properties of disk galaxies. It is con- 
venient to adopt power-law approximations for the initial 
virialized halo mass profile and angular momentum distri- 
butions, such that mh{£) ~ and /(£) - for small £. 
Figure 1 shows the location of various fiducial models in 
the plane of these power law indices m and n; the value 
m = 1 corresponds to the singular isothermal sphere, the 
value m = 4/3 corresponds to the Hernquist (1990) and to 
the Navarro, Frenk & White (1997) profiles, while the value 
m = 3/2 corresponds to the non-singular isothermal sphere 
with a constant-density core. These profiles span the range of 
dark-halo profiles suggested by theory, and plausibly consis- 
tent with observations. The shaded region is the allowed pa- 
rameter space for these models, constrained by the surface- 
density profile, rotation curve and disk-to-halo central mass 
ratio. The line ABC is the maximum angular-momentum in- 
dex consistent with a disk surface density profile in the cen- 
tral regions that declines with increasing radius, while the 
line DEF is the minimum angular-momentum index consis- 
tent with a finite value of the central circular velocity. The 
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1 4/3 3/2 

holo index m 



Figure 1. The power-law approximations for tlie initial virial- 
ized halo mass profile and angular momentum distributions axe 
mh{i) ~ i'^ and f{e) ~ ^" for small t The shaded region is 
the allowed parameter space for these models, constrained by 
the surface-density profile, rotation curve and disk-to-halo central 
mass ratio. The line ABC is the maximum angular-momentum in- 
dex consistent with a disk surface density profile in the central 
regions that declines with increasing radius, while the line DEF 
is the minimum angular-momentum index consistent with a fi- 
nite value of the central circular velocity. The line BF denotes 
the maximum values of n consistent with a non-zero central disk- 
to-halo mass ratio. The value m = 1 corresponds to the singular 
isothermal sphere, the value m = 4/3 corresponds to the Hern- 
quist (1990) and to the Navarro, Frenk & White (1997) profiles, 
while the value m = 3/2 corresponds to the non-singular isother- 
mal sphere with a constant-density core. 



line BF denotes the mciximum values of n consistent with a 
non-zero central disk-to-halo mass ratio. 

This power-law approximation has an initial virialized 
halo density profile at small radius as ph,ini{R) ^ -R^-™ 
(seen by solution of (22) for the form of g). Solving equations 
(21) - (23) above gives the corresponding profile after adiar 
batic infall. In the central region, where the disk dominates 
the gravitational potential (i.e. cof{£) ^ ^mh{£)), the halo 
density profile is ph{R) ~ R^^" ~ . Note that in the region 
where the dark halo dominates the gravitational potential 
(i.e. cof{t) ^irihit)), the halo density profile is essentially 
unaffected by the disk, as expected. 

For the case n = rn, the central halo density profile is 
unaffected since the disk mass density profile and the initial 
virialized halo density profile have the same dependence on 
i. The viable models within the shaded region have n < m, 
so that the final halo density profile in the central, disk- 
dominated region should be steeper than its initial virial- 
ized profile in this region, not surprisingly. Thus the final 
halo profile for these models ranges from ph ~ to 
Ph ~ R~^. It is interesting to note that ph ~ J?""'^^, the 
outcome of an initial virialized halo with a constant den- 
sity core responding to the settling of a disk with angular 
momentum index n = 4/3, corresponds to the de-projected 
de Vaucouleurs central density profile. Thus provided the 
virialized halo does not have a declining density profile with 
decreasing radius, which is unphysical, the final dark halo 
cannot have a constant density core, at legist in the very cen- 
tral region where the disk dominates, but should be cuspy. 



2.2 The Singular Isothermal Sphere 

The singular isothermal sphere provides a virialized halo 
density profile that is the most tractable analytically, and we 
can obtain some important scaling relations without having 
to specify the angular momentum distribution /(€); aspects 
of the analysis of this profile should hold in general, and 
provide insight. 

The final disk and halo mass profiles are given by solu- 
tion of: 



X 

rrih 
ma 



1 + 



(29) 

(30) 
(31) 



where < ^ < 1 and cq and ^, the parameters describing the 
compactness of the collapsed disk and its collapse factor, are 
defined above in equations (20) and (13). 

The collapse factor, defined as the ratio of the pre- 
collapse radius r2oo to the cutoff disk radius rc aX £ = 1 
(to be distinguished from the final disk scale length) is 



F)(l + Co) = 



V2\ 



V2 A 



(32) 



where c/, defined in equation (3), is a measure of the shape 
of the angular momentum distribution, and small values of 
Cf mean steeply-rising angular momentum distributions. 

Typical angular momentum distributions for disks are 
shown in Figure 2 (a,b) and discussed below. Here one 
should just note that a value c/ ~ 0.5 is reasonable. With 
this, and A — 0.06 and F — 0.1, we obtain ^ = 5.9 and 
Co = 0.66. The collapse factor as defined here is then about 
a factor of 9. For fixed angular momentum distribution c/, 
the collapse factor depends not only on 1/A but also on 
compactness co oc F/X. This is consistent with the results of 
previous studies of the collapse factor in two extreme cases 
(Jones & Wyse 1983; Peebles 1993): if the final disk is so 
self-gravitating that the value of F/X corresponds to Co ^ 1, 
then the collapse factor is cx F/X?; on the other hand if the 
final disk is sufficiently far from self-gravitating, with small 
F/X and co <C 1, the collapse factor is oc 1/A. The above 
collapse factor relationship is valid over the entire range of 
values for the parameters A and F (and is the most general 
relation derived to date). 

It should be noted that varying the normalized angular 
momentum distribution by varying c/ changes the derived 
collapse factor; this is investigated further below. 

Comparison between the sizes of observed disks and 
those predicted from such collapse calculations provides a 
constraint on the redshift at which the collapse happened 
(e.g. Mo, Mao & White 1998). Our disk cutoff radius for 
non-trivial f{t) may be expressed in terms of the initial con- 
ditions as : 

^=-10/7(.,)cHl-F+^f)' ^''^ 

where Zf is the 'formation' or assembly redshift, at which 
the halo is identified to have given mass and circular veloc- 
ity, with no mass infall after this epoch. One can see from 
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Figure 2. The angular momentum distribution function for the 
Mestel profile, f{t) = 1 — (1 — i)^/'^, is shown by the dashed line, 
(a) The angular momentum distribution function is }{b,tj = (1 + 
b)E-bf for models on line FED in Figure 1. The different curves 
correspond to 6 = 0, 1/4, 1/2, 3/4, 1. (b) The angular momentum 
distribution function is f {b, £) = (1 + lOfc)^*/^ - lOb^^z/is fo^ 
models on line BC in Figure 1. The different curves correspond 
to 6 = 0,0.1,0.3,0.6,1. 



this relation that both A and F/\ arc equally important; 
previous determinations considered the baryonic mass frac- 
tion fixed (Dalcanton, Spergel & Summers 1997; Mo, Mao 
& White 1998). 

The explicit inclusion of the parameter c/ allows us to 
take account of disk evolution, as gas is transformed to stars. 
Adopting a disk cutoff radius at three disk scale-lengths 
(e.g. van der Kruit 1987), and choosing specific values of the 
present-day stellar disk scale-length rd = 3.5kpc, A = 0.06, 
V2()() = 200 kms~^ and an Einstein-do-Sittor Universe with 
Hubble constant 0.5 < h < 1, the above relation gives the 
formation redshift 1.6 < Zf < 3.1 for F = 0.05 and c/ = 1/3; 
1.4 < Zf < 2.8 for F = 0.1 and Cf = 1/3; 0.9 < «/ < 2.0 
for F = 0.05 and c/ = 1/2; 0.7 < Zf < 1.7 for F = 0.1 and 
Cf = 1/2. 

Consistent with previous calculations, for fixed forma- 
tion redshift smaller values of A can lead to smaller disk size; 
we have here explicitly demonstrated that larger F/X can 
also lead to this result. Lower values of c/ also lead to higher 
redshift of formation. Thus for no viscous evolution i.e. there 
is no angular momentum redistribution during the evolution 
of the galactic disk and c/ has a time-independent value, the 
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formation redshift Zf determined by assuming fixed initial A 
and F and fixed disk size Tc is smaller than would be deter- 
mined if Cf could be decreased (to mimic say viscous evolu- 
tion). Lower values of c/ for given total angular momentum 
content imply a larger disk scale-length; the effect of viscous 
evolution is to re-arrange the disk material so as to increase 
the disk scale- length. Typical values of viscosity parameters 
lead to a factor of 1.5 increase in disk scale-length in a Hub- 
ble time. As we demonstrate below, a smaller value for c/ 
corresponds to larger bulge-to-disk ratio. This result is con- 
sistent with the results of Mo, Mao & White (1998): halo 
and disk formation redshift can be pushed to higher value 
when a bulge is included. These trends are general, and not 
tied to the specific model of the halo density profile. 

One should bear in mind that the old stars in the local 
thin disk of the Milky Way have ages of at least 10 Gyr, and 
may be as old as the oldest stars in the Galaxy (Edvardsson 
et al. 1993); the age distribution of stars at other locations 
of the Galactic disk is very poorly-determined, but it is clear 
that a non-trivial component of the thin disk was in place at 
early times (at redshift z > 2 for the cosmologies considered 
above). A common assumption in previous work is that the 
mass angular momentum distribution of the disk is that of 
a solid-body, rotating uniform density sphere, /(£) = 1 — 
(1 _£)3/2 (Mggtel 1963). For this distribution, c/ = 0.4, and 
one derives a low redshift of formation for a galaxy like the 
Milky Way (Mo, Mao & White 1998), which has difliculties 
with the observations. 



2.2.1 The Singular Isothermal Halo with Simple /(£) 

Analytic solutions to equations (29) - (31) can be obtained 
by assuming a simple monotonic increasing function /(6, 1) 
containing one parameter h with < & < 1. In order to avoid 
the situation where one obtains a trivial collapse factor due 
to a very small amount of disk material at very large radius, 
we restrict the shape of the mass angulax momentum dis- 
tribution function f[t) to avoid too shallow an asymptotic 
slope as / approaching 1 when I increases (see Figure 2). 
From Figure 1, the singular isothermal halo at point F re- 
quires f((.) ~ ^ for £ <C 1. A simple form consistent with this 
is f{t) = {l+b)£-M^ with < ^ < 1 and < 6 < 1. The an- 
gular momentum parameter c/ — (3 — &)/6. Figure 2a shows 
this angular momentum distribution function with difi'erent 
values of the parameter b, compared with the mass angular 
momentum distribution of a solid-body, rotating uniform 
density sphere, /(£) = 1 - (1 - £)^/^ (Mestel 1963). The 
normalized total angular momentum c/ corresponds to the 
area above each curve. 

Within the disk, where < £ < 1, the galactic disk 
surface density, circular velocity and the disk-to-dark mass 
ratio as a function of radius are: 

^ W\-2i>oU{z)F{L- Ffe 

^d[l) = -r, 

{1 + b- 2b£){l +C0 + cob - cob£f 

2£{1 + CO + cob) ' ^ ' 

Vc{£) = F2oo(l -F)[l + co(l + b - b£)], (35) 
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Ma{l) _ {l + b-M)(3-b)F 



6V2A(1 - F) 



where 

[l + c„(l + 6)]C(l 



F)R 



(36) 



(37) 



l + bco^{l- F)R ■ 

Prom equation (27)- (28), the circular velocity at radii be- 
yond the disk cutoff, but within the haJo, is: 



Vc{R) = 



V2oo{l-F) 



1 + ^1 + 



4F 



(1 - F)^R 



(38) 



(39) 



These results are plotted in Figure 3 (a,b,c) for A = 
0.06, F = 0.1, < £ < 1. The different curves correspond 
to 6 = 0, 1/4, 1/2, 3/4, 1. Larger values of b yield larger disk 
cut-off radii; note that the circular velocities for points be- 
yond the cut-off radius of a given model may bo obtained 
by forming the envelope of the values for the cut-off radius 
for larger values of 6. Varying the value of the parameter b 
changes the angular momentum distribution function simi- 
lar to the effects of viscous evolution; the ratio of disk mass 
to dark halo mass increases at small radius with increasing 
b, which can be interpreted as due to radial inflow of disk 
material. 

2.2.2 Halo Density Profile and Angular Momentum 
Combinations 

For initial virialized halo profiles other than the singular 
isothermal sphere, the collapse factor has a more general 
form: 

rc 

where rrihc is the solution of 1 = ^^rnhcg~^{rnhc) and is the 
mass fraction of the dark halo that is contained within the 
cut-off radius of the disk. The trend of the dependence of 
the collapse factor on the values of A and F remains the 
same as found above for the singular isothermal halo. 

The range of viable models represented by the shaded 
region in Figure 1 can be investigated by the appropriate 
virialized halo profile g{R) and angular momentum distri- 
bution function f{b, £) corresponding to the points E,B,D 
and C. The distribution of surface density, circular velocity 
and disk-to-halo mass ratio for these models, varying param- 
eter 6, are shown in Figures 4-7, by choosing fixed A = 0.06, 
F = 0.1 (and halo core size c = 4 if the halo has a core 
radius) . 

Model E: The results of a model corresponding to point 
E in Figure 1 are shown in Figure 4; this has a Hernquist 

halo profile g{R) = and f{b,£) = (1 + b)£ - b£'^ . 

The different curves correspond to 6 = 0, 1/4, 1/2, 3/4, 1. 

Model B: The results of a model corresponding to point 
B in Figure 1 are shown in Figure 5; this has a Hernquist 
halo profile g{R) = ^■^^^^ and f{b,€) = {1 + lQb)i'^'^ - 

IQbf'^^^^. This f{bj) is shown in Figure 2b. For I <^ 1, 
f{b = 0, f) ~ ^*/^; the 22/15 index in the second term is de- 
termined by the requirement that f(h,£) should be a mono- 
tonic increasing function of t, for all values of < & < 1. The 
different curves correspond to 6 = 0, 0.1, 0.3, 0.6, 1. 
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Figure 3. The model corresponding to point F in Figure 1. (a) 
the surface density, (b) circular velocity and (c) disk-to-halo mass 
ratio with A = 0.06 and F = 0.1. The virialized halo is a singular 
isothermal sphere. The angular momentum distribution function 

is ^) = (1 -I- b)l — b£^. The different curves correspond to 6 = 
0, 1/4, 1/2, 3/4, 1. Larger values of h yield larger disk cut-ofi' radii. 
The disk-to-halo mass ratio increases slightly when approaching 
the centre. 



Model D: The results of a model corresponding to 
point D in Figure 1 are shown in Figure 6; this has a 
non-singular isothermal halo with a constant density core, 
9iR) = 'T-ZTatf ' and f(b,£) = {l + b)£-b£\ Again the 
different curves correspond to fe = 0, 1/4, 1/2, 3/4, 1. 

Model C: The results of a model corresponding to point 
C in Figure 1 are shown in Figure 7; this has a non- 
singular isothermal halo with constant density core, g{R) = 
"'IZZZiT ' and fibj) = (1 + 106)£^/^ - 106£22/l^ Again 
the different curves correspond to 6 = 0, 0.1, 0.3, 0.6, 1. 
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Figure 4. The model corresponding to point E in Figure 1. (a) 
the surface density, (b) circular velocity and (c) disk-to-halo mass 
ratio with A = 0.06 and F = 0.1. The virialized halo is Hernquist 
halo with core size c = 4. The angular momentum distribution 
function is f{h, £) = {l + h)£ — bi'^ . The different curves correspond 
to 6 = 0,1/4,1/2,3/4,1. The disk-to-halo mass ratio increases 
significantly when approaching the centre. 

As can be seen from the figures, the different halo pro- 
files and angular momentum distributions produce disks 
with a variety of surface density profiles and rotation curves. 
As in Figure 3b, the circular velocities for points beyond the 
cut-off radius of a given model may be obtained by forming 
the envelope of the values for the cut-off radius for larger 
values of b. Thus if the disk is very compact, from equa- 
tions (27)- (28) or equation (38), we find that the circular 
velocity beyond the edge of the disk tends to decreases with 
radius. The location of the edge of the disk depends on both 
A and F (in addition to 6). Thus we have shown that rota- 
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Figure 5. The model corresponding to point B in Figure 1. (a) 
the surface density, (b) circular velocity and (c) disk-to-halo mass 
ratio with A = 0.06 and F = 0.1. The virialized halo is Hernquist 
halo with core size c = 4. The angular momentum distribution 
function is f{b, £) = {l + 10b)£^/^ - 106^22/15 -pj^g different curves 
correspond to 6 = 0,0.1,0.3,0.6,1. The disk-to-halo mass ratio 
increases but not significantly when approaching the centre. 

tion curves should show an imperfect disk-halo 'conspiracy' 
if the disk is too compact or too massive. This is consistent 
with observations (Casertano & van Gorkom 1991). 

Further, these figures demonstrate that with increas- 
ing b, the inner rotation curves become flat, and the tran- 
sition between disk-dominated and halo-dominated sections 
of the rotation curve becomes more and more smooth with 
increasing b. As discussed above, and illustrated in Figure 2, 
a higher value of b corresponds to a flatter specific angular 
momentum distribution function, and increasing b mimics 
the effects of viscous evolution in transporting angular mo- 
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Figure 6. The model corresponding to point D in Figure 1. (a) 
the surface density, (b) circular velocity and (c) disk-to-halo mass 
ratio with A = 0.06 and F = 0.1. The virialized halo is non- 
singular isothermal halo with constant density core size c = 4. The 
angular momentum distribution function is f(b, £) = (l+b)l—M'^. 
The different curves correspond to 6 = 0,1/4,1/2,3/4,1. The 
disk-to-halo mass ratio increases significantly when approaching 
the centre. 

mentum. This indicates that viscous evolution can help the 
creation of an apparent disk-halo 'conspiracy'. 

As can be seen from equations (21) - (26), for given 
Co (the compactness parameter defined in equation (20)), 
or F/\ ratio, within a given model of /(£, 6) and virialized 
dark halo profile g{R) , the normalized properties of the disks 
formed for 6 = are very similar. In particular the disk sur- 
face density profile, rotation curve and disk-to-halo mass 
ratio profile scale similarly with i. For example, in the case 
of the singular isothermal sphere (model F), equations (34) 
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Figure 7. The model corresponding to point C in Figure 1. 
(a) the surface density, (b) circular velocity and (c) disk-to- 
halo mass ratio with A = 0.06 and F = 0.1. The virial- 
ized halo is non-singular isothermal halo halo with core size 
c = 4. The angular momentum distribution function is /(fe, £) = 
(1 -I- 10fe)£*/^ - 106^22/15^ rj^Yic diff'crent curves correspond to 
b = 0,0.1,0.3,0.6,1. The disk-to-halo mass ratio increases sig- 
nificantly when approaching the centre. 

- (37) show that for 6 = the normalized disks are identical 
for given F/X ratio. Thus F/X must be an important factor 
in distinguishing one disk from another. The overall normal- 
ization of the surface density is oc F/X^ , so that again F/X 
and A enter separately and are both important. 

Note that here we are not insisting that the surface 
density profile of the gas disk so formed be exponential, 
unlike previous work (Mo, Mao & White 1998). Wc shall 
however appeal to viscous evolution tied to star formation 
to provide a stellar exponential disk. We now turn to this. 
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3 THE VISCOUS EVOLUTION AND STAR 
FORMATION 

In this paper we aim to link viscous evolution within disks 
and the Hubble sequence of disk galaxies. One of the mo- 
tivations for invoking viscous disks is that if the timescale 
of angular momentum transport via viscosity is similar to 
that of star formation, a stellar exponential disk is natu- 
rally produced independent of the initial gaseous disk sur- 
face density profile (Silk & Norman 1981; Lin & Pringle 
1987; Saio & Yoshii 1990; Firmani, Hernandez & Gallagher 
1996). Angular momentum transport and associated radial 
gas flows (both inwards and outwards) cam also, as shown 
above, provide a tight 'conspiracy' between disk and halo ro- 
tation curves, and, as demonstrated below, provide a higher 
phase space density in bulges as compared to disks. 

The star formation rate per unit area in a disk can be 
represented by a modified Schmidt law involving the dy- 
namical time and the gas density (Wyse 1986; Wyse & Silk 
1989). We shall use the form of the global star formation 
rate per unit area, E^, of Kennicutt (1998), based on his 
observations of the inner regions of nearby large disk galax- 



(40) 



where Ejas is total gas surface density, flgaa is the dynami- 
cal time at the edge of the gas disk, and the normalization 
constant, related to the efficiency of star formation, has the 
value a = 0.017 (Kennicutt 1998). Note that observations of 
the star-formation rates in the outer regions of disk galaxies 
suggest that it is actually volume density that should en- 
ter the Schmidt law, rather than surface density, and since 
many (if not all) gas disks flare in their outer regions, equa- 
tion (40) will over estimate the star formation rate (Ferguson 
et al. 1998). This is beyond the scope of the present model, 
but should be borne in mind and will be incorporated in our 
future work. 

For our models here, the edge of the initial gas disk is 
where i = 1, R = Rc and thus ^Igas = f2c- For halo formation 
redshift Zf, we can obtain the relationship between global 
star formation timescale and the galaxy initial conditions in 
the general form using equation (1), (25) and (39): 

= aQc = aVc{rc)/rc 

= WaH{zf)i{l-Ff{^mhc + cof, (41) 

where again irihc is the solution of 1 = ^'^mhcg~^{mhc)-i and 
is the fraction of the dark halo mass that is contained within 
the cut-off radius of the disk. The parameters ^ and Co are 
defined in equations (13) and (20); F is the baryonic mass 
fraction in the initial density perturbation. 

In the case of the singular isothermal halo, this relation 
has a simple form: 



(42) 



The gas consumption timescale will be longer than the char- 
acteristic star formation timescale due to the gas returned 

by stars during their evolution and death. For a standard 
stellar Initial Mass Function, tg ~ 2.5t, (e.g. Kennicutt et 
al. 1994). 
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This modified Schmidt law is based on observations of 
the inner regions of nearby large disk galaxies, and simple 
theoretical principles. Assuming it holds at all epochs allows 
one to estimate the properties of present-day disks from the 
initial conditions of earlier sections in this paper. Let us 
assume that the dark halo is fully virialized at redshift Zf. 
In keeping with the spirit of hierarchical clustering, let us 
allow for some star formation that could have taken place in 
the disk, from an earlier redshift Zi , and that the total mass 
of the system could increase until Zf (although to maintain 
the thin disk, this accretion and merging must be only of 
low mass, low density systems). So for any time t or redshift 
z between Zi and Zf, 



dMg 

dt 

dt 



= F 



dt 



dM» 
dt 



Mg 



(43) 

(44) 



where Mtot is defined in equation (1) and M, is mass locked 
up into stars. In an Einstein-de-Sitter Universe, H{z) = 
Ho{l + zfl"^ and H(z)t = 2/3. We have 



dt 

where A 



BMg 

^""^ and B 



(45) 



.... ^ - |a^(l - Ff{\ + co)^ and 

^ and Co are expressed in terms of c/, A and F through 
equations (13) and (20) in section 2.1. Hence, from equation 



(1), MtotOCj^y 

Identifying the dark halo to have fixed V200, indepen- 
dent of redshift, leads to A also being a constant, and thus 

the total mass grows as Adtot oc t. This differs from the stan- 
dard solution of infall onto a point, Mtot oc t^^'^ (Gunn & 
Gott 1972). 

The solution to the above equation is then 



Ma = 



A 



1 + B 



1 + B 



a) 



-(1 + B) 



(46) 



where ti corresponds to the redshift «» of the onset of star 
formation. Thus at the halo formation redshift Zf, the disk 

gas fraction is 



1 + B 



3(l+B)/2 



l + B 



(47) 



Since B depends on c/, A and F, i.e. B oc ^(1 — F + 
^y)^, the value of the constant B may be evaluated for 
A = 0.06 and various reasonable values of F and c/ as: 
for F = 0.1, Cf = 1/3, B = 0.30; for F = 0.05, c/ = 1/3, 
B = 0.23; for F = 0.1, c/ = 0.5, B = 0.59; and for F = 0.05, 
Cf = 0.5, B = 0.41. For fixed A and F, small values of B 
correspond to small values of c/, and hence flatter specific 
angular momentum distributions. 

Thus for 1 + Zf < 2(l-h^;i), fg{zf) x 1/{1+B), typically 
Z 2/3. 

We assume that after Zf, there is no further infall, and 
that the gas in the disk will be consumed with characteris- 
tic timescale tg{zf). Thus the gas fraction of a typical disk 
galaxy at the present time is 



fg = fg{zf) exp [-5t{zf) /tg{z})] , 



(48) 
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where dt{zf) is the time interval between the halo formation 
redshift Zf and present time z = 0. 

Assuming an Einstcin-dc-Sittcr universe, wc have 
dt(zf) = to{l - Ho/H{zf)) and toHo = 2/3. Thus a typ- 
ical value for the present gas fraction of disk galaxies is: 



In 



B[{l + Zf) 



3/2 



With the approximation fg{zf) = 1/(1 + B), we have 



-In/g -ln(l + B) 



(49) 



(50) 



For the Milky Way Galaxy, if we adopt 5 x IO^Mq 
for the atomic HI gas, and 1.3 x W^Mq for the molecu- 
lar H2 gas (Blitz 1996; Dame 1993), then with an estimate 
of 4 - 6 X 10^° Mq for the total baryonic mass of the Milky 
Way depending on stellar exponential scale length (Dehnen 
& Binney 1998), we obtain gas fraction, the total gas mass 
includes 24% helium mass, fg ~ 15% or oven higher do- 
pending on the mass model of the Galaxy. As mentioned 
earlier, for fixed A and F, small values of B correspond to 
small values of c/, and hence flatter specific angulax momen- 
tum distributions. For small values of the parameter B, say 
B ~ 0.3, we obtain Zf ~ 2.4, while for large values of B, 
B ~ 0.6, we obtain Zf ~ 1.2. The larger value of Zf is pre- 
ferred, given what we know of the age distribution of stars 
in the local thin disk (e.g. Edvardsson eto/.1993). The effect 
of viscous evolution is equivalent to choosing small c/, i.e. 
small values of B. So the inclusion of viscous evolution can 
give relatively higher halo formation redshift, which is con- 
sistent with the results from the constraint on the redshift of 
formation that we obtained from considerations of the size 
of the disk. 

It should be noted that for fixed halo 'formation' red- 
shift Zf, the star formation timescale tZ^ ~ j{l — F + 
j^j;)^- Again, both A and F/X are important. As we dis- 
cussed in the previous section, the structure of the normal- 
ized disk depends strongly on F/X while the overall nor- 
malization depends strongly on A for fixed F/X. As we show 
later, the bulge-to-disk ratio also depends on the importance 
of A and F/X. Thus many aspects of the Hubble sequence of 
disk galaxies - star formation timescale, disk gas fraction, 
bulge-to-disk ratio - depend on both A and F/X. 

The star formation timescale derived above is indepen- 
dent of V200 , which at first sight is surprising, given that the 
Hubble sequence of disk galaxies has been interpreted as a 
sequence of star formation timcscales, relative to collapse 
times (Sandage 1986), and observations show that the Hub- 
ble type of a disk galaxies is broadly correlated with the disk 
luminosity (Lake & Carlberg 1988; de Jong 1995). However, 
V200 is not an easily-observed quantity. 

The present-day luminosity of a disk in our model can 
be written as 



FMt^til - h) _ F[l - fgiB,Zf)]Vi 



7* 



10Gi/(«/)7, 



(51) 



with 7, the current value of the rnass-to-light ratio. Estima- 
tion of the predicted Tully-Fisher relation depends on what 
model parameter we use for the width of the HI line. Ob- 



viously if we identify this with V200 we have a reasonable 
relationship provided that the coefficient 7*/ is constant, 
i.e. _F[1 — fg{B,Zf)]/H{zf) — 7tf. This then requires that 
there should be a correlation between F and 2:/, in the sense 
that for large F, Zf is large. Then we can see from equation 
(42) that this leads to the star formation timescale being 
small, implying more efficient viscous evolution, and lead- 
ing to larger B/T ratio. Large zj may be correlated with 
small V200 in the context of hierarchical-clustering cosmolo- 
gies, and in that case, a short star formation time and large 
B /T ratio should be correlated with low luminosity, which is 
not consistent with observation. However, an interpretation 
that is compatible with observations is that high n-sigma 
fiuctuations for fixed ^200 can form high luminosity disks 
with large B/T ratio. 

The Tully-Fisher relationship is not a simple relation 
between luminosity and V200, but depends on where the cir- 
cular velocity Vc{R) is measured (Courteau 1997). The rela- 
tionship between V200 and Vc{R) obviously depends on the 
details of the halo density profile and the angular momen- 
tum distribution. From the rotation curves in fig. 3 to fig. 
7, we can see that it is appropriate to choose for our esti- 
mate of Vc the circular velocity at the cutoff radius of disk, 
adopted as 3 three scale lengths. From equation (35), we 
have Vc = 1^200 (l--F)(l + co), where co = c/F/^A(l-F), 
as defined in equation (20), is the compactness of disk. Now 
the predicted Tully-Fisher relation is 



L = 



F[l- }g{B,Zf)]V! 



10Gi/(«/)7*(l-Fco)3(l-F)3' 



(52) 



Requiring that the coefficient 7tf in the Tully-Fisher re- 
lation be constant, i.e. 7tf = F\l — fg{B,Zf)]/H{zf){l + 
co)^(l — F)^, then gives, from equation (42), oc co/(l -I- 
co). So for small co, the star formation timescale is large, 
which can cause less efficient viscous evolution. So less ef- 
ficient viscous evolution and small F/X will lead to small 
B/T ratio. Also small co is correlated with large Zf from 
the constancy of the coeflicicnt in the Tully-Fisher relation. 
Similarly large Zf may be correlated with small V200 in the 
context of hierarchical clustering cosmology. So small co and 
small V200 will lead to small Vc and lower disk luminosity. 
Thus this version of the predicted Tully-Fisher relation ap- 
pears fully compatible with the observations. However, one 
should boar in mind that we have adopted a fixed constant 
of proportionality a in the star formation law, and this may 
well vary with global potential well depth (White & Frenk 
1991) or local potential well depth (Silk & Wyse 1993). 

A further test of the model is the relation between disk 
scale and circular velocity, and its variation with redshift; 
observations indicate that Rd/Vc is smaller at high redshift, 
2 ~ 1 (Vogt et al. 1996; Simard et al. 1999). In our model 
there is little change in total mass between redsliifts of unity 
and the present, and so this evolution of disk size cannot be 
due to the halo mass growth, as had been proposed by Mao, 
Mo & White (1998). Instead in our model this is due to the 
different and changing scale lengths of gas and stars. In our 
model, due to early star formation, the gas fraction at Zf is 
about fg{zf) ~ 1/(1 -|- B), typically 2/3. It is natural that 
the gas component of the disk will have a larger scale length 
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than the stellar component. The scale length of the gas com- 
ponent can increase with time due to viscous evolution while 
the scale length of stellar component can also increase with 
time due to the non-linear local star formation law (Saio 
& Yoshi 1990). So the stellar scale length at high redshift, 
when the gas fraction is large, should be much smaller than 
the stellar scale length at present time, when gas fraction 
is lower. The study of detailed evolution of gas and stellar 
component is beyond the scope of this paper; but the pre- 
diction of stellar size evolution of our model is qualitatively 
consistent with the observations. 

The measured distribution of Rd/Vc for the local disk 
galaxy sample is approximately peaked at Rd/Vc ^ —1.5 
{Rd in kpc, Vc in kms~^) with a spread from -2 to -1 (Mao, 
Mo & White 1998; Courteau 1996). Our estimation of the 
disk size or scale length in Section 2 is valid for galaxies 
at the present time, when the gas fraction is small. Then 
assuming the disk cut off radius is three scale lengths, from 
equation (41), Rd/Vc = Rc/SVc = atf/3, and the present 
Rd/Vc is an indication of the galactic global star formation 
timescalc; further, from the constancy of 7tf, the coefficient 
of the TuUy-Fisher relation in equation (52), we have 

Rd/Vc = aU/3 = ^^TFJ. ],'^^° - (53) 

3 C0(1 - fg) 

Thus this ratio is also an indication of the disk compact- 
ness. Adopting the B-band mass-to-light ratio of our lo- 
cal disk 7* ~ 2.5Mq/Lq, and using the luminosity of our 
galaxy Li3 ~ 3 x 10^°L© and Vc ~ 220 kms~^ to esti- 
mate 7TF = Lb/Vc, we can obtain that \og{Rd/Vc) ^ 
-2.0-Hog( ^^^+_"»^) ). Obviously log{Rd/Vc) ~ -2 roughly 
corresponds to the predicted lower limit of the local sample, 
which is consistent with observations. The distribution of 
F from 0.05 to 0.2 and the distribution of A from 0.03 to 
0.12 will cause the value of the compactness parameter co to 
spread from 0.1 to 2 approximately. Adopting typical values 
F = 0.1, A = 0.06, Cf = 1/3, the peaJj will be located at 
\og{Rd/Vc) = —1.5 with Co = 0.44, which is again consistent 
with observations. The spread of Rd/Vc is simply caused by 
the spread in the value of compactness parameter cq. 



4 THE FORMATION OF BULGES AND THE 
HUBBLE SEQUENCE 

Galaxies classification schemes beised on morphology are 
the basic first step in understanding how galaxies form and 
evolve (van den Bergh 1998). The bulge-to-disk luminosity 
ratio is one of the three basic classification criterion for the 
Hubble sequence (Sandage 1961). However, the relation be- 
tween bulge-to-disk ratio and Hubble type has a fair amount 
of scatter, some of which must be related to the difficulty of 
a decomposition of the light profile into bulge and disk, and 
the bulge-to-disk ratio is dependent on the band-pass used 
to define the luminosity (de Jong 1996). The current obser- 
vational data show that bulges can be diverse and heteroge- 
neous (Wysc, Gilmorc & Franx 1997). Some share properties 
of disks and some are more similar to ellipticals. Models of 
bulge formation can be classified into several categories: the 
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bulge is formed from early collapse of low angular momen- 
tum gas, with short cooling time and efficient star formation 
(Eggen et al. 1962; Larson 1976 - who invoked viscosity to 
transport angular momentum away from the proto-bulge; 
van den Bosch 1998); the bulge is formed from merging of 
disk galaxies (Toomre & Toomre 1972; Kauffmann, White 
& Guiderdoni 1993); the bulge is formed from the disk by 
secular evolution process after bar instability (Combes et al. 
1990; Norman, Sellwood & Hasan 1996; Sellwood & Moore 
1999; Avila- Reese & Firmani 1999); the bulge is formed from 
early dynamical evolution of massive clumps formed in disk 
(Noguchi 1999). However it is speculated that large bulges, 
which tend to have dc Vaucouleur's law surface brightness 
profiles, share formation mechanisms with ellipticals, while 
smaller bulges, which tend to be better fit by an exponen- 
tial profile, are formed from their disks through bar dissolu- 
tion. It should be noted that the significantly higher phase 
space density of bulges when compared to disks suggests that 
gaseous inflow should play a part in the instability (Wyse 
1998), just as invoked in earlier sections of this paper, for 
other reasons. 

Here we only consider this latter case, the formation of 
small bulges from the disk. Early studies of disk instabilities 
showed that Toomre's local stability criterion Q > 1 is also 
sufficient for global stability to axisymmetric modes (Hohl 
1971; Binney & Tremaine 1987). It is known that the bar in- 
stability requires a similar condition (Hockney & Hohl 1969; 
Ostriker & Peebles 1973). Efstathiou, Lake & Negroponte 
(1982) used N-body techniques to study the global bar in- 
stability of a pure exponential disk embedded in a dark halo 
and proposed a simple instability criterion for the stellar 
disk, based on the disk-to-halo mass ratio. However, it has 
been argued recently that there is no such simple criterion 
for bar instability (Christodoulou, Shlosman & Tohline 1995; 
Sellwood & Moore 1999). Further, recent N-body simula- 
tions (Sellwood & Moore 1999; Norman, Sellwood & Hasan 
1996) show that every massive disk form a bar during the 
early stages of evolution, but later the bar is destroyed by 
the formation of a dense central object, once the mass of 
that central concentration reaches several percent of the to- 
tal disk mass. This can be understood in terms of the linear 
mode analysis work and nonlinear processes of swing am- 
plifier and feed back loops (Goldreich & Lynden-Bell 1965; 
Julian & Toomre 1966). Toomre (1981) argued that the bar- 
instability can be inhibited in two ways: one is that a large 
dark halo mass fraction can reduce the gain of the swing 
amplifier, while the other is that feedback through the cen- 
tre can be shut off by a inner Lindblad resonance (ILR). 
The dense central object can destroy the bar via the sec- 
ond of these mechanisms. However, most of these studies 
assume that the dark halo has a constant density core. On 
the contrary, we have shown that dark halo profile cannot 
have a constant density core after adiabatic infall, if one 
starts from physical initial conditions. It will be interesting 
to study bar formation under different disk-halo profiles in 
addition to the well-studied harmonic core. 

We will for simplicity adopt the simple criterion of Ef- 
stathiou, Lake & Negroponte (1982), interpreted to deter- 
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mine the size of a bar-unstable region, with the radial extent 
of the bar, ri defined by 



Md{n)/Mh{n) > 13, 



(54) 



with the value of the parameter /3 chosen to fit observations. 

As wc have argued in section 2, F/X is the important quan- 
tity determining the overall normalization of the disk sur- 
face density. From N-body simulations (Warren et al. 1992; 
Barnes & Efstathiou 1987; Cole & Lacey 1996; Steinmetz 
& Bartelmann 1995), the distribution of A can be fit by a 
log-normal distribution: 



P{X)d\ ■■ 



exp 



ln2(A/Ao) dA 



2ax 



A 



(55) 



where Ao = 0.06 and a\ = 0.5 (this result is fairly indepen- 
dent of the slope of the power spectrum of density fluctua- 
tions). 

What of the possible range of the baryonic fraction F? 
Some previous studies suggested that _F ~ A as an explana- 
tion of the disk-halo 'conspiracy' (Fall & Efstathiou 1980; 
Jones & Wyse 1983; Ryden & Gunn 1987; Hernandez & 
Gilmore 1998). Some interpretations of the observed Tully- 
Fisher relationship suggest that indeed F is not invariant 
(McGaugh & de Blok 1998). Here we shall assume that the 
distribution of F, similar to A, is log-normal, centreed at 
Fo = 0.1 with ffF = 0.05. So F is mainly within the range 
0.05 ~ 0.2. 

We generated a Monte Carlo sample of disks, with fixed 
halo V20Q and halo formation redshift Zf, but with values of 
F and A following the above distributions. Then for given 
virialized dark halo profiles g{R) and angular momentum 
distribution f{b,i) in Figure 1, we can calculate the star 
formation timescale from equation (41). The parameter b 
represents the efficiency of viscous evolution; on the assump- 
tion that the viscous timescale is equal to the star formation 
timescale, we can use a simple linear correlation between the 
value of b and the star formation timescale. Then the value 
of parameter b can be obtained. From equations (21), (22), 
(23), (26), (50) and (54), one can calculate the bulge-to-disk 
ratio and final disk gas fraction. 

Thus we can plot B/T ratio versus disk gas fraction, 
for the whole Monte Carlo sample, confining the parameter 
space of A and F to give small bulges with B/T < 0.5. 
Larger values of B/T would correspond to such an unstable 
disk that the exercise is invalid. Disks that are too stable, 
having low F/\ or large A, will evolve little and probably 
end up as low surface brightness systems. 

Figure 8a shows model F, representing the singular 
isothermal halo. Assuming the typical values A = 0.06 and 
F = 0.1, then the choice of /3 ~ 0.8 is required to allow 
the existence of bulges. This is just the value of /? used in 
the global bar instability criterion (Efstathiou, Lake & Ne- 
groponte 1982). For this model the disk-to-halo mass ratio 
varies little with radius (as shown in Figure 3a), so the B/T 
ratio is strongly dependent on \/F. Only a small range of 
values of A/F is allowed, so as to not over-produce either 
bulge-less disks or completely unstable disks. 

The relation between B/T ratio and gas fraction for 
this model is given in Figure 8b; there a general trend in 
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Figure 8. The model corresponding to point F in Figure 1. (a) 
The sample of disk galaxies with different A and F. (b) The rela- 
tion between B/T ratio and gas fraction for these sample galaxies 
at current age. The different curves in this plot correspond to dif- 
ferent values of (1 — F)\/F A/F, which sets the overall trend. 



the observed sense, with large scatter. The different curves 
in this plot correspond to different values of A/F, which sets 
the overall trend. 

Figure 9 (a,b) shows the equivalent plots of model E, 
representing a Hernquist profile halo model, which is prob- 
ably a more realistic case. Here the the disk-to-halo mass 
ratio varies strongly with radius when approaching the cen- 
tre, and the allowed parameter space for values of A and of F 
that allowing the formation of bulges can be large. Overly- 
unstable disks are denoted by Eisterisk symbols in Figure 9a, 
and low surface-brightness disks are denoted by cross sym- 
bols. The relation between B/T ratio and gas fraction is 
similar to that for the isothermal halo. 



4.1 Constraints from the Milky Way 

The Milky Way bulge is reasonably well-fit by an exponen- 
tial profile, with a scale-length approximately one-tenth that 
of the disk (Kent et al. 1991). The morphology of the bulge 
is consistent with some triaxiality (Blitz & Spergel 1991; 
Binney et al. 1997). Perhaps the Milky Way is a system in 
which the bulge has formed from the disk, through a bar 
instability? Observations show no evidence for a significant 
young or even intermediate-age stellar population in the field 
population of the Galactic Bulge (Feltzing & Gilmore 1999), 
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Figure 9. The model corresponding to point E in Figure 1. 
(a) The sample of disk galaxies with different A and F. Overly- 
unstable disks are denoted by asterisk symbols. Bulgeless disks 
or low surface-brightness disks are denoted by cross symbols, (b) 
The relation between B/T ratio and gas fraction for these sample 
galaxies at current age. The different curves in this plot corre- 
spond to different values of (1 — F)\/F A/F, which sets the 
overall trend. 



despite their being ongoing star formation in the inner disk. 
This imphes that if the bulge were formed from the disk 

through bar dissolution, only one such episode is allowed, 
and this should have happened at high redshift. In the eon- 
text of the present model, the lower star formation rates 
and longer viscosity timescales of later times act to stabilize 
the system. However, it remains to be seen if the observed 
relative frequencies of bars, bulges and central mass concen- 
tration is consistent with the models of bar dissolution. 



5 SUMMARY 

In the context of hierarchical clustering cosmology, the dark 
halo of a disk galaxy can be formed by quiescently merging 
small sub-halos into the primary dark halo, or by smoothly 
accreting matter into the dark halo. We derive the generic 
solution to the adiabatic infall model of disk galaxy for- 
mation pioneered by many authors (Mestel 1963; Fall & 
Efstathiou 1980; Gunn 1982; Faber 1982; Jones & Wyse 
1983; Ryden & Gunn 1987; Dalcantori, Summer & Spergel 
1997; Mo, Mao & White 1998; Hernandez & Gilmore 1998). 
Through exploring the allowed parameter space of dark halo 



Formation of Hubble Sequence 13 

profile and angular momentum distribution function, we 
show that the central halo density profile should be cuspy, 
with the power law index ranging from —0.75 to —2, in the 
central regions where the disk mass dominates. 

Using a modified Schmidt law of global star formation 
rate, we derive a simple scaling relationship between the disk 
gas fraction and the assembly redshift. We explicitly allow a 
distribution in the values of the baryonic mass fraction, F, in 
addition of the distribution in values of the spin parameter 
A. These two are found to play different role in determining 
the structural properties of the final disk, the star formation 
properties and bulge-to-disk ratio. 

We mimic viscous evolution of disks by varying the spe- 
cific angular momentum distribution of the disk, to redis- 
tribute angular momentum as a function of time. We de- 
rive a consistent picture of the formation of galaxies like the 
Milky Way, with old stars in the disk. Under the assumption 
that the viscous evolution timescale is equal to the star for- 
mation timescale, we can further combine the A and F with 
the efficiency of angular momentum redistribution caused 
by viscosity. Assuming that small bulges are formed from 
their disks through bar dissolution, we can use the global 
bar instability condition to obtain bulge-to-total ratio, and 
explore the dependence on F, X and the viscous evolution 
efficiency. 

The inclusion of viscous evolution has the merits of ad- 
dressing several important issues: the conspiracy between 
disk and halo, the formation of the exponential profile of 
stellar disk, the high phase space density of bulges. We have 
presented an analytic treatment, to illustrate these points 
and identify areas of particular need for more work. 
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